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Abstract 

c7^ ' We analyse deterministic diffusion in a simple, one-dimensional setting consisting of a 

04 , family of four parameter dependent, chaotic maps defined over the real line. When iterated 

under these maps, a probability density function spreads out and one can define a diffusion 
coefficient. We look at how the diffusion coefficient varies across the family of maps and 
under parameter variation. Using a technique by which Taylor-Green-Kubo formulae are 
V,^ ' evaluated in terms of generalised Takagi functions, we derive exact, fully analytical expres- 

pH , sions for the diffusion coefficients. Typically, for simple maps these quantities are fractal 

functions of control parameters. However, our family of four maps exhibits both fractal 
and linear behavior. We explain these different structures by looking at the topology of 
the Markov partitions and the ergodic properties of the maps. 

PACS numbers: 05.45.Ac, 05.45. Df, 05.60.Cd 
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^ ' 1 Introduction 

en 

• ' Chaotic diffusion lies at the heart of the interaction between transport theory and dynamical 

Li^ , systems [T]-[5]. The study of this interaction was started in the early 1980's [3HS] and since then 

f— ^ ' defines an active area of research producing interesting results for both mathematicians and 

physicists [7l|8]. By analysing transport theory in the context of dynamical systems, one can 
take into account the deterministic equations of motion and therefore does not need to rely on 
any stochastic based analysis. 

In this paper, we evaluate the diffusion coefficients for a set of piecewise linear, one- 
dimensional chaotic dynamical systems. There are a few methods for analytically evaluating 
the diffusion coefficient of a dynamical system; there is the cycle expansion technique described 
in [S], there is an interesting method based on the zeros of the systems dynamical zeta func- 
tion [7], there is the ^ first passage method^ |10] from statistical physics which is based on the 
escape rate of the system and involves computing the eigenmodes of the probability density 
and there is the ^twisted eigenstate method^ described in [3]. We are interested in how the 
diffusion coefficient varies under parameter variation. The diffusion coefficient is often a very 
complicated fractal function of the parameter in simple one dimensional maps, this result first 
being reported in [TTl . Unfortunately, analytically the above methods can often only be worked 
out for specific values of control parameters, or are otherwise rather difficult to handle. We 
will thus use a method based on the Taylor-Green-Kubo formula [IHS] and certain fractal 'gen- 
eralised Takagi functions^ from |10| . in order to highlight this methods power at dealing with 
the complicated structure that one typically finds in the diffusion coefficient. Furthermore, to 
this end we choose our parameter in such a way as to simplify the derivation of the diffusion 
coefficients relative to the work in jTO. Our main results are the analytic expressions for the 






parameter dependent diffusion coefneients of these maps, an explanation for the mixture of hn- 
earity and fractahty in the diffusion coefficients and the discovery that the diffusion coefficients 
are actually very stable in certain parameter ranges. 

In section[2]we introduce the family of maps that are under scrutiny before briefly describing 
the method we use to obtain the parameter dependent diffusion coefficients in section [31 In 
section 2] we discuss and define functional recursive relations in terms of infinite sums for the 
interesting fractal Takagi functions that we meet. In section [5] we analyse the structure of the 
diffusion coefficients and explain the features that we observe, most importantly we explain why 
we get a mixture of fractality and linearity. We also employ a method based on the Markov 
partitions of the maps to pinpoint exactly where the local extrema will be in these fractal 
diffusion coefficients [TUHT^ . In addition we explore the fact that in certain parameter regions, 
the diffusion coefficients are very stable to dramatic changes in the microscopic dynamics of the 
maps, this property being interestingly opposed to the fact that the diffusion coefficients are 
extremely sensitive to parameter variation in other areas of the parameter space. Section [6] is 
a conclusion. 

2 The family of maps 

In this section the family of maps and how they are constructed will be introduced. 

For h > let Mfi{x) : [0,1] -^> R be a parameter dependent variant of the well known 
Bernoulli shift map. The parameter lifts the first branch, and lowers the second branch i.e. 

M'^(^) = [ 2x-l-h i<x<l • W 

In order to create an extended system for diffusion, we define Mh{x) : M — > M by periodically 
copying equation ([T]), with a lift of degree one such that 

Mh{x + n) ^ Mh{x) + n, n e Z. (2) 

The map described above and its diffusion coefficient were first studied in [13j using a 
different method. We call it the lifted Bernoulli shift map. This process of copying a map 
with a lift of degree one is a common way to create a diffusive map [IHS] . The use of the lift 
parameter h ensures that the invariant probability density function (p.d.f) remains a constant 
function throughout the entire parameter range. This helps simplify the derivation of the 
diffusion coefficient as an invariant p.d.f is an essential ingredient in the Taylor-Green-Kubo 
formula that is used to derive the diffusion coefficient [UITO] . 

The remaining members of the family are created by changing the sign of the gradient in 
equation IQ. Let Whix) : [0, f ] ^ R 

...,.( -2x + h + l 0<x<^ 

W'^(^) = [ -2x + 2-h i<:r<l ' ^^^ 

which we call the lifted negative Bernoulli shift map. Let Vh{x) : [0, 1] ^ M 

^'^(^)- 2x-l-h h<x<l ' (4) 
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which we call the lifted V map. Let Aii{x) : [0, f] 



^>^(^)-[ ^2x + 2^h i<.<f ' ^^) 

which we call the lifted tent map. Again we apply the lift of degree one condition of equation 
(HI) to equations Q, (01) and ([5]) to create spatially extended systems defined over the real line. 
See Figure [T] for an illustration. 
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Figure 1: The family of maps. In this figure, a section of eaeh of the four maps under serutiny is illustrated at 
a parameter value of h = 0.5. In (a) the lifted Bernoulli shift, so-called because the Bernoulli shift is recovered 
when h = on the unit interval. In (b) the lifted negative Bernoulli shift, so-called because a version of the 
Bernoulli shift with a negative gradient is recovered when h = 0. In (c) the lifted V map, so-called because a V 
map is recovered when h = 0. In (d) the lifted tent map, so-called because a tent map is recovered when h = 0. 



3 Deriving the diffusion coefficient 

The diffusion coefficient D is given by Einstein's formula for diffusion in one dimension 



D{h) = lim 



2n 



(6) 



where we are evaluating the diffusion coefficient as a function of the parameter h. The angular 
brackets (...) represent an average taken over the invariant p.d.f, p*{x) 



(...) :- / dxp*{x). 
Jo 



(7) 



Equation ([5]) can be transformed into the Taylor-Green-Kubo formula [T] 



D{h) = lim [Y, {io{x)ik{x))] - \ {io{xf) (8) 

\k=0 / 



where Vk{x) 



Vk{x) := Xk+i -Xk (9) 

gives the displacement of a point x at the fc*'* iteration Xk- It can in fact be shown that we 
need only consider the integer value of the displacement of a point x [H] , i.e. we can replace 
Vk{x) with the simpler function Vk{x) : M — > Z, 

Vk{x) := \xk+i\ - [xk\ , (10) 

which for the lifted Bernoulli shift map Mh{x) is 



[/ij <Xk < 



l-h 



vk{x)^{ i'- 74 ' (11) 

-[h\ ^<ik<l 

where [...J and [...] are the floor function and ceiling function respectively and x is a; modulo 
1. The function /i : M — > [0, 1] is defined as 

h=l 1 ^^^^> . (12) 

[_ h mod 1 otherwise 

Equation ()12p is simply a corrective function that ensures equation ()11[) is correct at the points 
of discontinuity. We can also use the fact that p* (x) — 1 in our model to simplify equation ([5]) 
further to 

D{h) = lim ( / vo{x)Yvk{x) dx ] - ^ I VQ{xfdx . (13) 



/ vo{x)'Yvk{x) dxj - - 



Evaluating equation (I13p will give us our diffusion coefficient. The second integral is simple 
enough, however the first integral is more involved. We will now explain the method we use for 
this. 

We first define a cumulative ^jump function^ J^.;{x) : [0, 1] — ;■ R as 

n 

J^Ax):^Y.''k{x) (14) 

which gives the integer displacement of a point x after n iterations. The subscript M tells us 
we are considering the jump function for the lifted Bernoulli shift map Mfi{x). We substitute 
equation ([Ti| into equation (IT5|) to obtain 

DM{h)= lim ( / vq{x)J'1;j{x) dx] - - / vl{x) dx. (15) 

n->oo \J^ J 2 Jq 

In order to extract the information we need from p3|) . we define Tm{x) : [0, 1] ^- M which 
integrates over the jump function, as in the first integral of (J15l) . 



Tm{x) :- lim Tliix) = lim / J^Mdy- (16) 

71— >00 71— >-00 J f^ 

Equation (fTH]) defines the 'generalised Takagi functions' discussed above and in [10]. Due to the 
chaotic nature of the maps, and in particular, the sensitive dependence on initial conditions, the 



jump function behaves very erratically for large n. This is reflected in Tm{x) which becomes 
fractal in the limit as n — )■ oo. We call these functions ^generalised Takagi functions^ because 
for the lifted Bernoulli shift map with parameter h = 1, one obtains the function first studied 
by T. Takagi in 1903 [15]. This function is also an example of a 'De Rham' function |16) . 
T. Takagi was interested in this function from an analytical perspective as a function that is 
both continuous and non-diffcrentiable. 

In order to control the chaotic nature of the jump function, we derive a functional recursive 
relation [T^. For the lifted Bernoulli shift map 



Jl,{x) ^ v^{x) + J^^-' [aIu{x)) 



(17) 



where Mh{x) is equation ([T]) taken modulo 1. In turn, we define a functional recursive relation 
for the generalised Takagi functions, by substituting equation (J17p into equation pi 



where tM{x) : [0, 1] 



Tm{x) = lim Tlf{x) 

n— J-oo 

= lim [tM{x) 
M is defined as 



It;,-! (M,(a 



(18) 



tnix) 



vo{v)dy 



xvq{x) + c. 



(19) 



The constant of integration c is evaluated by requiring the continuity of the Takagi function 
and that r(0) = T{1) = 0. We also use the symmetry of the Takagi functions to simplify the 
evaluation. For the lifted Bernoulli shift map we obtain the recursion relation 



Tm{x) 



\Tm hx + h] + [h\x- ^TM{h) 



0<x < 



l-h 



/^] a; + %i - ^TM{h) 1^ < x < ^ 



\Tm {2x-h^+^^- \h-\ X - \TM{h) 
\Tm [2x^l-h\- \Ji\ X + \Ji\ - \TM{h) 



2 

^ <x<l 



(20) 



where we have taken the limit n — > oo. We now have the ingredients that we need to derive the 
parameter dependent diffusion coefficient. We firstly apply equation (|lip to the Taylor-Green- 
Kubo formula to obtain 



DmQi) 



lim 

n— >-oo 



vq{x)JIi{x) dx 





l-h 



2 7o 



Vq[x) dx 



= lim / [h\Jli{x) dx+ \h]J2Ax) dx 



[h\Jli{x) dx - 




\h]Jlj{x) dx (21) 



\hr 



\hr 



lh\' 



Evaluating the integrals, simplifying using equation (I20|) and gathering relevant terms we obtain 



DMih) = {[h\-\h^)(^[h\^hlh\-TM{h^) + \h](^\h]+h-l^ 



(22) 



-([h\'{i-h) + \hrh], 



which after some wranghng can be rewritten as 



l2 



DM{h)=^-^+l^jil^2\h]) + TM{h) . (23) 

Equation (|23p is an analytic expression for the parameter dependent diffusion coefficient for 
the lifted Bernoulli shift map. The first two terms in equation p3| form a piecewise linear 
function that is asymptotically equal to -^ for large h. This term tells us how the diffusion 
coefficient grows for large h. Interestingly, our shifted map belongs to a different universality 
class compared to the maps studied in |10[I17] where the gradient was varied as a parameter. 
There it was conjectured that the parameter dependent diffusion coefficient for this class of 
maps would grow quadratically with the height parameter h with a coefficient of i. 

The Tm ( h j term tells us about the fine scale structure. This fine scale structure is periodic 

modulo 1 as it is a function of h, see figure [5] for an illustration. We also observe a second 
region of asymptotic behaviour where the diffusion coefficient tends to ^ as h ^- 0, i.e. it 
behaves linearly and reproduces a simple random walk result for diffusion. This result agrees 
with the findings in [IHIIIZ], it tells us that the higher order correlations of our system are 
negligible as ft. — > 0. This change from one type of asymptotic behaviour to another in a system 
was denoted as a crossover in deterministic diffusion [171 . 

Although the calculations have been presented for the lifted Bernoulli shift, the method is 
the same for the other three maps. If we turn our attention to the lifted negative Bernoulli 
shift map Wh {x) the Takagi function is given by 



Tw{x) 



^Tw(-2x + h) + \h]x + ^Tw{h) 0<x<~^ 



-iTn- (-2a; + h + l) + [h\ x + '^ {\h'] - [ftj) + ^Tw{h) ^<x<\ 



-^Tw ( -2a; + 2 - ftj - [/ij a; + [/ij + | ([ft] - [ftj) + \Tw{h) 5 < a; < 1 - | 
'\Tw (-2x + 3 - ft) - [ftl a; + [ft] + \Tw{h) 1 - | < 2; < 1 

(24) 
and the corresponding expression for the parameter dependent diffusion coefficient is, 

Dwih) - ^ + ^ {\hf - [ft]') +Tw{h). (25) 



The first two parts of equation ([251) form the same piecewise linear function found in equation 
(j23p so as ft ^- 00 we observe the same asymptotic behaviour found in the lifted Bernoulli shift 
map, see figure [H 

The parameter dependent Takagi functions for the lifted tent map and lifted V map are 
different in character to the two Bernoulli shift maps discussed above. In order to emphasise 
this difference we restrict the parameter to ft e [0, 1]. For the lifted tent map the Takagi function 
is 

iTA(2a; + ft) - \Tj~,{h) < .x < i^ 

, a; + irA(2a; + ft-l) + ^-iTA(ft) ^ < 2; < i 

^^^"^^"^ -irA(-2:E + 2-ft) + | + irA(I-ft) i<a;<I-^ ' ^^^^ 

-X - iTA(-2a; + 3 - ft) + 1 + iTA(I -ft) I-|<a;<i 

Whilst the Takagi function for the lifted V map looks like 



x-\Tv{-2x + h) + \Tv{h) 0<a;<| 

-iry(-2.T + I + ft) + | + iTy(ft) h<x<\ 

-a;+iTv(2a;-ft) + | + i-iTy(I-ft) l<x<^ 

\Tv{2x-l-h)-\Tv{l-h) li^<x<l 



(27) 
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Figure 2: Large scale structure and asymptotic behaviour. In this figure, the diffusion coefficients for the hfted 
BernoulU shift (a) and the lifted negative Bernoulli shift (b) are illustrated. Also included in both is /(h) = h'^/2 
to show how the function grows for large h. Note the periodicity of the fine scale structure. Inset in both is an 
illustration of the asymptotic behaviour as /i ^ O.The dashed lines arc ^ to show the different behaviour in the 
two maps. 

The thing to notice here is the inclusion of the T{1 ~ h) terms. Usually we would use the 
symmetry of the Takagi functions to equate these terms with T{h), but the Takagi functions 
for the lifted V map and lifted tent map do not have this symmetry. We will discuss the 
consequences of this and the diffusion coefficients for these two maps in sections |4] and El 

4 Evaluating the Takagi functions 

In this section we will show how to evaluate the Takagi functions and see what some of them 
look like in figure [31 We will start with the Takagi functions for the lifted Bernoulli shift map 
(equation (PU)) ). We can repeatedly apply the recursion relation in equation P^ and obtain a 
sum 



Tm{x) 



tM{x) + -Tm{Mu{x)) 



-TMih) 



= tM{x) + \tM{Mn{x)) + -Tm{MI{x)) 

OO ^ 

fc=0 



-TM{h) - -TM{h) 



(28) 



where 



tmix) 



\h\x 
%i + \K\ X 

[h\ - [h\ X 



<x < 



i-/t 



l-h 



< X < 



\<x<i4 

^^<x<\ 



(29) 



We can remove the Tu^h) term in equation 
(HI 



to obtain an infinite sum in terms of equation 
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Figure 3: The Takagi functions. In this figure, the Takagi functions are shown for the four maps at a parameter 
value of h = 1. In (a); the hfted BernouUi shift, which you may recognise as the famous Takagi function. In 
(b); the negative BernoulH shift. In (c); the hfted V map. In (d); the hfted tent map. Note the self similarity 
and 'fractal' structure. Note also the asymmetry in (c) and (d) compared to (a) and (b), this is due to the 
asymmetry in the evolution of the p.d.f for these maps. Also portrayed here is that Tv{x) = Ta(1 — x), this 
result will be explained. 



oo 



fe=0 



E ^tM (m,^(x)) - ^tm (h) - Me 4*m f^''('^) ) - ^^e^) 



fc=0 



\k=0 



°° 1 / 1 \ 



k=0 



(30) 



The Takagi functions for the hfted negative Bernoulh shift map can be evaluated in a similar 
way. We use the same recursive definition to obtain 






k=0 



k=a 



E^ twix) + -Tw{h) 



fc=0 



(31) 



where 



\h']x < a; < I 
twfx)-=<' lm-[h\)+[h\x |<x<i 

^ [h\ + I (pi] -[h\)~lh\x i < X < 1 - I ■ 

\h]-\h']x 1 - I < a; < 1 
We can again remove Tw{h) to obtain an infinite sum in terms of equation p2p . 

r^(x) = f;(^^) [t^{wi;{x))+\tw(wi;{h))y (33) 

Things are not so simple when we try to evaluate the Tagagi functions for the lifted V map 
and the lifted Tent map. The fact that we have both a positive and a negative gradient in 
Vh{x) and Kh{x) makes it harder to simplify the recursion relation into one single sum. If we 
take the Takagi functions for the lifted V map as an example and let 

m--[-\ ?|:^! , (34) 

nix) -I ^^('^^ 0<x<i 

then the sum for the Takagi function that we obtain is 

Ty{x) = ty{x) + \g{x) + y^ ^^^\y^^ (tviV^ix)) + \9{Vl:{x))] (36) 

where 



2^ ' ^-^ 2^ \v\h\ii 2- 

fc=i 



X <a; < I 
tv{x):={ ^_| tl!;!:L ■ 0</i<l. (37) 

2 2 — ■^ *^ 2 

^^ <X<\ 

So this difficulty is the first consequence of the asymmetry in these maps. We can however 
simplify the Takagi functions for the lifted V map and the lifted tent map by setting h = \ and 
hence evaluate them. See figure [3] for an illustration of some of these interesting functions. 

5 The structure of the diffusion coefficients 

In this section the parameter dependent diffusion coefficients for the four maps will be shown, 
and their structure explained. 

5.1 The lifted Bernoulli shift map 

Figure m (a) gives the diffusion coefficient for the lifted Bernoulli shift map. The two striking 
features are the fractal region when h is between zero and a half, and the linear plateau when 
h is between a half and one. These regions will be explained in turn. 

5.1.1 The fractal region 

Firstly, the term ^fractaP has no strict mathematical definition so we use the term loosely. In 
particular we use it to refer to the fact that the diffusion coefficient exhibits non-trivial fine 
scale structure, and regions of scaling and self similarity; for a discussion of this see [TH]. 

The topological instability of the map under parameter variation is refiected in the fractal 
structure of the diffusion coefficient. So in order to understand the fractality, we need to 
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Figure 4: The diffusion coefficients. In this figure, the parameter dependent diffusion coefficients are illus- 
trated. In (a) the lifted Bernoulli shift. In (b) the Ufted negative Bernoulli shift. In (c) the lifted V map. In 
(d) the lifted tent map. 

understand the topological instability. To this end, we take equation ([l} modulo 1, and analyse 
the behaviour of the Markov partitions of the interval map Mh{x) : [0, 1] — > [0, 1] 



h 0<x<^ 



Mh{x) 



2x 
2x + h~l 

2x ~h 5 < x < 
2x-l-h 



2 — ^ "^ 2 

1 ^ ™ ^ 1+h 



(38) 



The structure of the Markov partitions of (j38]) varies wildly under parameter variation. The 
method we employ to understand the Markov partitions involves iterating the critical point 
X = ^, see [TDHH]- The set of iterates of this point, along with the set of points symmetric 
about x ~ ^, will form a set of Markov partition points for the map. Hence we call the orbit of 
X = ■^ a 'generating orhiV . Furthermore, if the generating orbit is finite for a particular value 
of h, we obtain a finite Markov partition. We can then use the finite Markov partition to tell 
us about the diffusive properties of the map and hence the structure of the diffusion coefficient. 
For this purpose the following proposition is crucial. 
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Proposition 1 The set of values of the parameter h which give a finite Markov partition are 
dense in the parameter space. 

Proof: 

We show that when h is rational, the generating orbit is finite. This is achieved by showing 
that the denominator of h fixes the number of possible iterates of the generating orbit. Let 
h — ^ where a,b G N and a < b 

Ma {0.5) ^ 1-h 
b — a 



(39) 

Clearly, b — a £ {0, 1, 2, ..., 6 — 1}. If Mn. (^^) is then evaluated, there are four possibilities due 
to the four branches of ([55]) . It is not hard to show that for all four possibilities 

Mi[^-^]=T^ ce{0,l,2,...,b-l}. (40) 



So for n = 1 and n = 2 



6' 
with c G {0, 1, 2, ..., b — 1}. Now we assume that 



6' 



MUQ.5)^l (41) 



MT (0.5) = 7 (42) 

'' b 

for all m < n, and the result follows by induction that 

M2(0.5) = 7, ce {0,1,2,...,6-1} VneN. (43) 

The result in (|43l) puts a limit on the size of the subset of values that the orbit of x — 0.5 
can hit at a given rational value of h. This size being equal to |{0, ^ ..., ^^}| = b. Hence the 
orbit must be periodic or pre-periodic, and the Markov partition of the map must have a finite 
number of partition points when h is rational. 

q.e.d. 

The second important result is that the finite Markov partitions correspond to the local 
minima and maxima of the diffusion coefficient |10fll2j . If we extend our view back to the full 
maps, we see that if the generating orbit is periodic, i.e. 

M;^{0.5) = 0.5, n e N, (44) 

then this corresponds to a relatively high rate of diffusion for the parameter value, which is 
reflected in the diffusion coefficient as a local maximum. In contrast, if the generating orbit is 
pre-periodic, this corresponds to a relatively low rate of diffusion for the parameter value, which 
is reflected in the diffusion coefficient as a local minimum. So given that we have a dense set 
of local maxima and minima, we observe a fractal diffusion coefficient. Furthermore, equation 
(|44p furnishes us with the means of pinpointing the local maxima of the diffusion coefficient. 
In [lOj this technique gave an approximation of where the local extrema were, but this model 
allows to to find them precisely. Each n gives us a set of simple linear equations for the variable 
h, the solutions of which are the local maxima in the diffusion coefficient. In addition, the 
smaller values of n give the most striking local maxima. One can apply a similar technique to 
locate the local minima of the graph, see figure [S] However, it's not a case of finding solutions 
to one simple equation like (j44[) . Rather, there are many ways to define pre-periodic orbits as 
opposed to defining periodic orbits. In addition, the local minima do not adhere to such a strict 
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ordering that the local maxima do. This is due to the fact that there are two components to a 
pre-periodic orbit, namely the transient length and the periodic orbit length. In summary, for 
h between zero and one half, there exists a dense set of points which are either local maxima 
or local minima. Hence a fractal structure is observed. 



D(h) 




Figure 5: Pinpointing the local extrema. In this figure, some of the local extrema have been highlighted. At 
the points highlighted by squares we see the orbit of 0.5 is iterated to infinity resulting in a local maximum. At 
the points highlighted by diamonds we see the orbit of 0.5 is in a closed loop resulting in a local minimum. 



5.1.2 The linear region 

The second feature of the diffusion coefficient is the linear region where 0.5 < h < 1. Not only 
is it striking because it very abruptly changes from fractal to linear, it is also counter intuitive 
if we apply a simple random walk approximation to the map. That is, we can obtain a first 
order approximation of the diffusion coefficient by looking at the measure of the escape region 
of the map (the area where a point can move from one unit interval to the next) [101117] . This 
region clearly increases linearly with the parameter, so based on a first order approximation one 
would expect to see a general increase in the diffusion coefficient as the parameter increases. 
However, what is observed defies this. We explain this feature by noting the non-ergodicity of 
the map. 

When the parameter h reaches one half, a fixed point is born in the modulo 1 map. As 
the parameter increases further, the fixed point bifurcates and the two resulting fixed points 
split the phase space up into two invariant sets, breaking the ergodicity of the map, see figure 
[6l Consequently, the invariant density p*{x), can be interpreted as the sum of two invariant 
densities pl{x) and P2ix) and the diffusion coefficient can be evaluated as 
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Figure 6: Nonergodicity. In this figure, the nonergodicity of the the hfted BernouUi shift map at h = 0.75 is 
shown in (a), and in the hfted negative BernouUi shift at /i = 0.25 in (b). For simphcity the dynamics have 
been reduced to the maps modulo 1. One can see that the black areas get mapped into themselves as do the 
grey areas. This splits up the phase space breaking ergodicity. 



1 f^ 

n->oo 2n Jq 



= lim -— 

n— )-oo 2n 



pl{x){Xn ~ Xo)^dx + / P2{x){Xn - Xo)'^dx 

Jo 



(45) 



We can then use the Taylor-Green-Kubo formula and derive two separate generalised Takagi 
functions, in order to evaluate the diffusion coefficient as 



DMih) = {l-h) + ih--) 

_ 1 
2' 



(46) 



5.2 The negative Bernoulli shift 

Figure |3Kb) shows the diffusion coefficient for the lifted negative Bernoulli shift map Wh{x). 
Firstly we note how radically different the structure of the diffusion coefficient is from the 
lifted Bernoulli shift map. In addition, as /i — ?• the diffusion coefficient does not go to | by 
reproducing the random walk solution as in the lifted Bernoulli shift map. However, as /i ^> oo 
the diffusion coefficient goes to -^ as in the lifted Bernoulli shift. We also observe both a 
linear and a fractal region. For < h < ^ the diffusion coefficient is a simple linear function 
equal to h. Again, the explanation for this is that the map is non-ergodic in this parameter 
range, see figure [H The phase space is split up into two invariant regions, one of which does 
not contribute to diffusion as it is a trapping region, the other of which grows linearly with 
h so we observe a linear, increasing diffusion coefficient. When ■^ < h < 1, the map becomes 
topologically unstable under parameter variation. Similar to the lifted Bernoulli shift map, this 
instability is reflected in the behaviour of the Markov partitions of (|4]) taken modulo 1. We have 
that the finite Markov partitions are dense in the parameter space and that these finite Markov 
partitions correspond to the local maxima and minima of the diffusion coefficient. Hence we 
observe a fractal diffusion coefficient. 
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5.3 The lifted V map 

Mercifully, we need not turn to infinite sums like equation p6p when we evaluate the diffusion 
coefficient for the lifted V map. The diffusion coefficient for the V-map is given by 

Dv{h) = \ + \ (Tvih) + TvH - h)) ,0 < h < I. (47) 

Instead of evaluating the relevant Takagi function, equation (|27p . numerically we can use the 
helpful property that for h less than one half 

Tv{h) = -lTvil-h) + f^ + ^Tv{h) 

= -Tv{l -h) + h. (48) 

Using equation pS)) in equation P7|. the diffusion coefficient can be evaluated to 

Dv{h) = ^ + ]^{-Tv(l-h) + h + Tv{l-h)) 

= h. (49) 

Furthermore, for h greater than one half 

Tv{h) = -h+]^Tv{h) + ^ + ]^-\Tv{\-h) 

= -Tv{l-h)-h+l. (50) 

Using equation ((50)) in equation (|47| , the diffusion coefficient can again be evaluated 
Dv{h) = ^ + ^{-Tv(l-h)-h + l + Tv(l-h)) 



(51) 



See figure [41(c) for an illustration. 



5.4 The Hfted tent map 

For the lifted tent map Ah{x), we can not perform the same trick with the Takagi functions, 
equation (pS)) . that we did with the lifted V map in subsection 15.31 However we still do not 
need to resort to numerical computations to see what the diffusion coefficient looks like as we 
can note that 

A,.(x) = ~Vhi-x). (52) 

Equation ([5^ is important because it serves as a topological conjugacy of the form /(x) = —x. 
By using the Taylor-Green-Kubo formula equation p^ it was shown in [19] that the diffusion 
coefficient is preserved under topological conjugacy, hence we observe an identical diffusion 
coefficient in the two maps which can be seen in figure SI An alternative way is to start from 
Einstein's formula equation ([S]), as is demonstrated by proving the following proposition: 

Proposition 2 The diffusion coefficients for the lifted tent map is identical to that of the lifted 
V map. 
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Proof: 

We note that the diffusion coefBcient for the hfted V map is given by 

Dv{h) = hm ^ ^, (53) 

n— >oo Zn 

and that the diffusion coefficient for the tent map (Z?a(/i)) can be given by 

{(A-i-xo) - {^xo)f) 
DA{h) = hm ^ '-, (54) 

n— >oo Zn 

where the p.d.f is in the interval, [—1,1] to make it symmetric about a; = 0. By substituting 
dUD into dMD 



DA{h) = hm 



(AJJ(-xo) - {^xo)f 



n—nx 2n 

2 



= hm 



i^Vf^^ixo) - i^xo)y 



n— >oo 2tI 

({V^{xo)-xof 
= hm ^ 

n— >oo 2n 

= Dv{h). (55) 

we arrive at the desired relationship that D\{h) — Dv(h). 

q.e.d. 

From what we have seen in the lifted Bernoulli shift map and the lifted negative Bernoulli 
shift map, we would expect to find non-ergodicity in the lifted tent and lifted V maps. Fur- 
thermore we would expect to find it across the entire parameter range given the linear diffusion 
coefficients. However, there is no obvious non-ergodicity, and although a proof that the maps 
are ergodic across the entire parameter range remains elusive, we can check for ergodicity at 
individual values of the parameter by checking the reducibility of the transition matrices [20j . 
So we can confirm ergodicity for some values of h. Furthermore, the behaviour under parameter 
variation of the Markov partitions of the maps Q and (O taken modulo 1 indicate the presence 
of topological instability in the parameter space. So, although all the ingredients for a fractal 
diffusion coefficient are present, we observe a linear one. 

In order to understand the linearity of the diffusion coefficient of the hfted V map (and 
hence the hfted tent map also) , we first note the similarity of the linear regions in the diffusion 
coefficients of the two lifted Bernoulli shift maps, see figured This presents the question, why 
do these maps have the same diffusion coefficients when they have such different microscopic 
dynamics? 

We will explain the linearity of the diffusion coefficients of the lifted tent and lifted V map 
by showing why they have the same diffusion coefficients as the two Bernoulli shift maps in 
the relevant parameter ranges. We will take these ranges in turn starting with 0.5 < h < 1. 
The diffusion coefficient for the lifted Bernoulli shift map is given by equation (|23p which for 
h e [0.5, 1] simplifies to 

DM{h) = ^+TM{h) (56) 



2 



which can be rewritten as 



h .. f f^ tM^ ^^ ^ tMih) \ 




DMih) = - + Jim I n ^^1 + -^ (57) 
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See Appendix (jX| for a full derivation of equation ([57)) which tells us how the diffusion coefficient 
converges as n — >■ oo. Now, keeping equation (j57|) in mind, we turn our attention to the lifted 
V map. The diffusion coefficient for the lifted V map with h G [0, f] is given by 

Dvih) = I + hm (]-T--\h) + ir{J-i(l - h)) . (58) 



From equation (|27p we can derive the useful recursion relation 



2 ^ ' 2 



r^(/i) = ty(h) + -T^-\h) - -r^-\\ - /i) (59) 



which if we continue to apply leads to 

n "1 

^vC^) = E ^i^e^) - E ^^r'(i - '^)- (60) 

k= 

Using equation ([60)1 in equation (|58| 



2fc ^ ^ ^ /l^ 2*^ 

fc=0 fc=l 



'• , ,:„ i'-^ h.(h) l' 



^-w-i+i™ Ei#+5ri'"<'-'"-i^^^'"""<'-''> ■ <•"' 



Now, we note that 



^fc=0 \ k=l 



n-l 



r^i(f-M = ^^r^i-'=(f-/.), n^oo. (62) 

fe=l 

so these terms cancel each other out in the limit and we are left with 

Dv{h)=^^+tv{h). (63) 



If we let n ^^ cxD in equation [S7] then 

tMJh) 
2" 
and hence we are left with the expression 



2 



n -^ cx) (64) 



DM{h) = ^+tM{h). (65) 

We can see from equation ([^^ and equation (1571) that tM{h) = tv{h) when /i G [0.5, 1] and hence 
the two diffusion coefficients are equal despite A'Ih{x) being non-ergodic and Vh{x) ergodic in 
this parameter range. However the two equations (|5T|) and ([571) tell us that DM(h) and Dv{h) 
converge at different rates, see figure ([7]) for an illustration of this phenomenon. We also note 
that the diffusion coefficient is only dependent on ij\/(/i) or tv{h) and that these functions are 
only dependent on the branch of the map in [0.5, 1]. We interpret this phenomenon physically 
as the diffusion undergoing a dominating branch process, i.e. the diffusion coefficient is only 
dependent on the contribution of one branch of the map in the limit as n — >■ cxd. Hence we see 
identical diffusion coefficients despite the different microscopic dynamics. 

We find a similar situation for h G [0,0.5], where the lifted negative Bernoulli shift map 
and the lifted V map have the same diffusion coefficient. They have the branch in [0,0.5] in 
common in this parameter range and it is this which creates the dominating branch process. 
Hence we also observe identical diffusion coefficients between these two maps in this parameter 
range. 
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20 40 60 80 

n 

- Negative Bernoulli shift V-map| 



Figure 7: The time dependent diffusion coefficients. In this figure we see how the diffusion coefficient converges 
for certain parameter values. In (a) h = 0.2 for the lifted negative Bernoulli shift map and the lifted V map, in 
(b) h = 0.7 for the lifted Bernoulli shift map and the V map. We see that the diffusion coefficients tend to the 
same value but at different rates, indicating that the difference in the microscopic dynamics docs not play a role 
in the limit. Rather we observe a dominating branch process where the common branch of the map determines 
the diffusion coefficient. 

5.5 Stability of D(h) in the non-ergodic regions 

In this subsection we look at the two hfted Bernoulh shift maps in the non-ergodic regions. We 
have aheady seen that changing the gradient of one branch of the map (resulting in the lifted V 
map or lifted tent map) has no effect on the diffusion coefficient in these regions, even though 
the microscopic dynamics are affected greatly. We will explore this phenomenon further. 

For h G [0,0.5] the lifted negative Bernoulli shift map has a linear diffusion coefficient. 
By changing the gradient and chopping up the second branch of the map we obtain a map 
Wh{x): [0,1]^R 



-2x 



Wh{x) = 



2x- I 
2x-2 



0<a;< i 

1 <a;<l-| 



(66) 



The Takagi function for this map with h G [0, 1] is 



Twi^) 



-\T^{-2x + h)+x+\T^{h) 

ir^(2a;-l) + | 
^T^{2x^l) + l-x 



0<a;< I 

2 
< X < I - 

_ h 
2 



I 
2 
I 



h 

2 

<X <l 



(67) 



and the diffusion coefficient can be evaluated as 



D^{h) ^\ + \ (T^ih) + T^{1 -h)),0<h<l. 



Again for h G [0, 0.5] we have that 



T^{h) = -T^il - h) + h. (69) 

which implies that for h G [0,0.5] the diffusion coefficient is equal to h. 

We can play a similar game with the lifted Bernoulli shift map which has a linear diffusion 
coefficient for h G [0.5, f]. We can change the gradient and chop up the first branch of the map 
to create Mhix) : [0, 1] ^ M 



Mh{x) 



-2a; + I 0<a;<i^ 
-2a; + 2 ^ <x <^ 
2x-l-h i<a;<I 



(70) 
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-\T^„{-2x + l) 




Q<x<^ 


-iT^^(-2a; + l)+x-i^ 






\T^{2x-h)-x+^-^-\T^,{l- 


-h) 


lT^^{2x-l-h)~\T^,{l-h) 




1+'' < rr < 1 



The Takagi function for this map with /i G [0,1] is 



T-l-r^—} 2-^My ^■^ ~^ ^> ~^ ■^ 2 2 —•■'-• ^ 2 (7^\ 



and the diffusion coefficient for this map is 

DmW = \ + \ (^M W +TM(l-h)),Q<h< 1. (72) 

We see that for h E [0.5, 1] we have 

T^j{h)^^T,-,{l-h)-h+l (73) 

which imphes that for h G [0.5, 1] the diffusion coefficient is equal to 0.5. 

So we have again seen that the diffusion coefficients for these maps are very stable in the 
relevant parameter ranges, i.e. the ranges where the maps are non-ergodic. As long as the 
diffusion coefficient is given by 

D{h)='^ + lTih) + ^Til-h), (74) 

we can manipulate one branch of the Bernoulli-shift or the negative-Bernoulli-shift maps and 
the diffusion coefficient will remain unaffected in the non-ergodic parameter ranges, despite the 
fact that the non-ergodicity may be broken. 

6 Conclusion 

We have derived exact analytical expressions for the parameter dependent diffusion coefficients 
of four one dimensional maps. This was achieved by using Taylor-Green-Kubo formulae and 
generalised Takagi functions. Under parameter variation we have observed both fractal and 
linear behaviour in the diffusion coefficients. The fractality was explained in terms of the topo- 
logical instability of the maps under parameter variation and this was understood by analysing 
the Markov partitions of the map. The linearity was explained in terms of the non-ergodicity 
of the maps in certain parameter ranges, this non-ergodicity splits the phase space up into 
two ergodic components each with their own diffusion coefficient. These individual diffusion 
coefficients compliment each other to create a linear diffusion coefficient. We also observed 
linear diffusion coefficients despite all the hallmarks of fractality being present like topological 
instability under parameter variation [3l ll0ffT2] and ergodicity. In this case we found that in the 
relevant parameter range, the ergodic maps have a set of branches in common with the non- 
ergodic maps, these common branches dominate the diffusion process in the long time limit and 
hence we observe identical diffusion coefficients. In addition, when the parameter causes these 
maps to be non-ergodic, we found that the diffusion coefficients of these maps are so stable that 
we can drastically alter the microscopic dynamics without affecting the diffusion coefficient. 
This finding serves as a counter example to the previously held belief that if your system was 
topologically unstable and ergodic, you would observe a fractal diffusion coefficient. It is also 
a counter example to the belief that if you have a linear diffusion coefficient then you have a 
topologically stable system. 

Future work will involve finding out under exactly what conditions we can manipulate 
the microscopic dynamics of these maps and still observe the same diffusion coefficient. We 
would also like to learn whether there exist any other systems which display this dominating 
branch phenomenon. We could potentially apply the techniques used here to more realistic, 
higher dimensional systems and see if we still obtain similar results. Also of interest is the 
consequences of introducing a bias into the system generating a current. It would be worthwile 
to study whether analogous phenomena exist for this other transport property and whether 
they can be revealed by similar techniques. 
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A Takgagi function recursion relation for Mh{x) 

We can modify the Takagi function recursion relation for the lifted Bernoulli shift map in order 
to understand how the diffusion coefficient converges. We restrict the parameter h to [0, 1] and 
see from equation pO| that 

TUh) = tM{h) + \Tlr\A'h{h))~\Tl-\h) 

= tM{h) + ^tMiMnih)) + -T-r\Ml{h)) \nr\h) - ]-T-i\h) (75) 



where we have not taken the limit n — )■ cxd. If we continue to apply equation ()20|) we arrive at 
the recursive definition 

^mW - E 2^^M {^^(^)) E ^^rr'(M' (^6) 

k=0 k=l 

In order to simplify equation (I76p further we define 



rH-.^Y.^nr'ih)- (77) 

fc=i 

We can write equation (j77p recursively as 

Substituting equations (1771) and ([75)1 into ([75)1 we obtain 



T{n)^-T;;:,-\h) + -T{n-l). (78) 



1 

Tnih) = E^^M(A^,t(/i))-r(n). 

k=0 
" 1 1 1 

= T.^iM{^<(''))-2^M\h)--Tin-l). (79) 

fc=0 

Then subbing equation ([75)) back into equation ([75]) 



n /n— 1 _. n — 1 _. \ _. 

fe=0 \fe=0 fc=l / 

= E ^*^^ (^^''e^)) - 2 E ^*A/ (Mtih)) + -r{n - 1) - -r(n - 1). (80) 

fc=0 \k=0 I 

We then arrive at the expression 

' -tM{Mim (81) 



nAh) = (e ^^m (a4^(M)) + ^' 

\fe=0 / 
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which is a useful expression for the Takagi functions as it only involves terms which contain 
tAiix). For 0.5 < /i < 1 we can simplify equation (j8ip further by using the fact that h is a fixed 
point of the modulo 1 map, Mfi{x) 



\fc=o / 



(82) 



so our expression for the diffusion coefhcient in the range 0.5 < /i < 1 is 



Dnih) = -+\imTlY\h) 



2 n—^oo 

S+ii-ffE^V^^V (83) 

2 n^ooM^-^ 2*^+1 / 2" / 
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